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Using a density functional method, we investigate the prop- 
erties of liquid ^He droplets doped with atoms (Ne and Xe) 
and molecules (SFe and HCN). We consider the case of 
droplets having a quantized vortex pinned to the dopant. A 
liquid drop formula is proposed that accurately describes the 
total energy of the complex and allows one to extrapolate the 
density functional results to large A'^. For a given impurity, 
we find that the formation of a dopant +vortex-|-*Hejv com- 
plex is energetically favored below a critical size Na- Our 
result support the possibility to observe quantized vortices in 
helium droplets by means of spectroscopic techniques. 



Since the first observation of the 1/3 vibrational band 
of SFg dissolved in ''He droplets 0, the infrared spec- 
troscopy of molecules inside or attached to helium has 
attracted a wide interest (see, for instance, Refs. |^-^] 
and references therein). A major motivation for these 
efforts is that cold helium droplets offer the possibility of 
resolving rotational spectra of rather complex molecules 
and may constitute 'the ultimate spectroscopic matrix' 
1^ to create and study novel species. This unique fea- 
ture of helium droplets originates from their quantum 
nature: not only they are fluid at zero temperature, due 
to the large zero point motion, but also exhibit a cru- 
cial superfluid behavior. The superfluid character of ^He 
droplets is interesting also from a fundamental viewpoint. 
In fact, the observation of superfluid effects in finite-sized 
quantum systems has to do with important concepts, like 
order parameter, Bose-Einstein condensation, and phase 
coherence, which were originally introduced for uniform 
systems and which are now widely used in different con- 
texts. 

In the case of liquid helium, Grebenev et al. recently 
showed that only a rather small amount of "^He atoms is 
needed to develop a superfluid droplet, also confirming 
theoretical predictions |^. In that experiment, the evi- 
dence for superfiuidity is the appearence of a sharp ro- 
tational spectrum of an OCS molecule in ■^He-^'He mixed 
drops, when the number of *He atoms surrounding the 
dopant is larger than about 60. In the same spirit, ex- 
periments have been made to observe critical velocities 



Q (i.e., the occurrence of a Landau criterion for super- 
fluidity), and a reduction of the moment of inertia (see 
pO| and references therein). In contradistinction, detect- 
ing quantized vortices in droplets still remains an open 
question. It is worth stressing that all these investiga- 
tions have many analogies with the current activity on 
Bose-Einstein condensation in trapped gases, where new 
results are now available about critical velocities [0 , mo- 
ments of inertia ||l2|,^ and vortices |lj] . 

In this work we address the problem of quantized vor- 
tices. One first observes that a vortex line in a pure 
droplet is expected to be difficult to produce and stabi- 
lize, since it implies a significant increase of energy com- 
pared to a vortex-free droplet. In order to circumvent 
this limitation, we explore the possibility of pinning the 
vortex line to a dopant atom or molecule. If the dopant 
is deeply bound inside the droplet, it might stabilize the 
vortex for a time long enough to permit its observation. 
A second advantage is that the dopant could make the 
detection feasible via spectroscopic techniques. 

Our purpose is to determine the energy and density 
profile of a impurity-|-vortex-|-*HeAr complex, for droplets 
up to N= 1000, using a finite-range density functional. 
We then subtract to its energy that of the same droplet 
without vortex and/or impurity and show that the differ- 
ence fits very well to a liquid drop formula, which allows 
one to safely extrapolate to larger droplets. The density 
functional method consists in minimizing the total energy 
of the system at zero temperature written as a functional 
of the He density. We use the Orsay-Paris functional , 
which is based on an effective non-local interaction with 
a few parameters fixed to reproduce known properties of 
bulk liquid He. This functional has been shown to accu- 
rately reproduce the static properties of pure and doped 
He clusters |Q and has also been used to describe a quan- 
tized vortex line in bulk liquid helium [ p7| . In the latter 
case, the vortex is included with the Feynman-Onsager 
(OF) ansatz for the velocity field. This implies a singular 
vorticity and hence the vanishing of the density on the 
vortex axis. At T= this approximation is a reasonable 
starting point, since it enormously reduces the compu- 
tational cost. Recent calculations have shown that 
the density profile and energy of the vortex line given by 
the OF approximation are reasonably close to the ones 
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obtained by assuming non-singular vorticity. Finally, the 
actual temperature of the droplets, below 0.4 K is low 
enough for neglecting thermal contributions to the fluid 
motion. Thus, corrections beyond OF and at T 7^ are 
expected not to change the main results of the present 
work. 

The minimization of the energy is performed in axial 
symmetry by mapping the density on a grid of points, 
putting the vortex line along the z-axis and the dopant 
in the center, at r = 0. The numerical code used to 
calculate the density profile and energy is the same used 
in [ p^ . The potentials for rare gas impurities have been 
taken from [|9| , that of the spherically averaged SFg from 
|§, and that of HCN from 

We first consider pure droplets with and without vor- 
tex. In Fig. 0we show density profiles, at z = 0, obtained 
for different N. For large droplets the shape approaches 
that of a rectilinear vortex in the uniform liquid [0 ; the 
core radius is of the order of 1 2 A and the density 
oscillates as a consequence of the He-He interaction. In 
Fig. H we plot the energy associated to the vortex flow, 
defined as 



AEy{N) EE Ev{N) ~ E{N) 



(1) 



where i?v and E are the energies of droplets with and 
without vortex, respectively. The solid line represents 
the results obtained with a liquid drop formula of the 
kind: 

AEy{N) = aN^/^ + (iN^I^ logiV + 7iV~i/3 , (2) 



with the parameters a= 2.868 K, /3= 1.445 K, and 7= 
0.313 K extracted from a fit to the density functional 
calculations. This formula works well. The reason can 
be easily understood by means of a hollow-core model for 
the vortex, having core radius a, in a droplet of radius 
R and constant density po- By integrating the kinetic 
energy of the vortex flow in the limit i? 3> a, one gets 
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Writing R— r^N^/^ one recovers the A^-dependence as in 

The next step is the inclusion of a dopant atom or 
molecule. As an example, in Fig. || we show the He den- 
sity distribution for a drop of iV = 500 with HCN hosted 
in the vortex core. Both the axis of the linear molecule 
and that of the vortex are taken along z. The density 
is very inhomogeneous near the dopant, due to the com- 
plexity of the HCN-Hc interaction. The energetics of the 
system can be conveniently analysed by introducing the 
following energies: 



AE^iN) 
Sx{N) 



Ex+v{N)-Ex{N) , 
Ex{N)-E{N) , 
Ex+viN)-E{N) , 



(4) 
(5) 
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where the subscripts X and V refer to drops doped with 
impurity X and/or vortex line. 

The energy A£^y is the one associated with the vortex 
flow in the doped cluster. In Fig. ^ it is compared with 
the vortex energy in pure droplets, AEy. The difference 



SxiN) = AE^ - A£:v < 



(7) 



is almost independent of N, apart from the smallest 
droplets. The reason is that this difference has to do 
with the 'geometrical extension' of the dopant, i.e., the 
'hole' made by the dopant in the vortex flow, as well as 
with the distortion of the density near the dopant caused 
by the pinning of the vortex core. Both effects are local- 
ized near the dopant and, thus, they are expected to give 
a shift in energy which becomes A'^-independent for large 
droplets. 

The quantity Sx{N) in Eq. is the solvation energy 
of the dopant in a vortex- free droplet. The results ob- 
tained for Ne, Xe, HCN and SFg are shown in Fig. ^ As 
already discussed in |jl^,^ , the solvation energy becomes 
almost iV-independcnt for N larger than a few hundreds. 
The value at = 1000 can be safely taken to represent 
the solvation energy in the bulk, Sx{oo) ~ 5x(1000). For 
our analysis, we have chosen impurities having binding 
energies on a wide range. 

The key quantity in the present study is the solvation 
energy of the dopant -I- vortex complex given by S'x+v(A^) 
in Eq. (^). The results are shown in Fig. ||. From the 
deflnitions (0)-(0) one can also write 

Sx+viN) = ExiN) + AE^iN) - E{N) 

= Sx{N) + AEv{N) + Sx{N) . (8) 

In Fig. ^ we compare ^x-i-v with the sum Sx + AEy; 
the difference is Sx- The simple picture which emerges 
from this analysis is that the solvation energy of the 
dopant -|-vortex complex is just the sum of the solvation 
energy of the dopant with no vortex and the extra en- 
ergy of a vortex in a pure droplet, apart from a small 
shift which depends on the dopant. Deviations from this 
rule are signiflcant only for small droplets, having ra- 
dius of the order of the size of the dopant. Our numer- 
ical results provide a quantitative basis for this picture 
and yield typical estimates of 5x- It is worth noticing 
that, by rearranging the terms in (|7|), this quantity can 
be written as the difference between the solvation ener- 
gies of the dopant in a droplet with an without vortex, 
6x = [(-E-x+v ^ Ey) — {Ex — E)], and can hence be inter- 
preted as the binding energy of the dopant to the vortex 

il; 

Since the solvation energy Sx is negative and almost 
constant for N > 300 while the vortex energy AEy al- 
ways increases, the dopant -f-vortex complex has a solva- 
tion energy which changes sign at some Na-. This means 
that for N < N^, the dopant -t- vortex complex is ener- 
getically favored. In the case of Ne, as one can see in 
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Fig. H, Nci ^ 380. This number is rather small as com- 
pared to the typical droplet size in current experiments, 
and is a consequence of the weak binding of Ne. Dopants 
with stronger binding have larger Nc^- An estimate for 
HCN, Xe, and SFg can be easily obtained by means of 
the liquid drop formula. One has to insert expression (|^) 
in ^ and use the large- iV values of Sx and Sx- The 5'x 
values turn out to be —310 K, —320 K and —622 K, and 
the Sx values are 5.0 K, 4.4 K, and 7.7 K, for Xe, HCN, 
and SFg, respectively. These numbers yield TVcr ~ 7600 
for Xe, ~ 8100 for HCN, and ~ 40000 for SFg. 

In conclusion, the analysis of the energetics of doped 
helium droplets has allowed us to disclose a possible 
mechanism to create and stabilize vortex lines. A 
dopant -|-vortex complex could be formed by picking up 
the impurity, assuming that the collision imparts suffi- 
cient angular momentum. The vortex line is expected to 
appear attached to the dopant, since the binding energy 
i5x is negative. The formation of the complex is energet- 
ically favored below a critical N which is well within the 
range of droplet sizes met in current experiments if the 
dopant has a large solvation energy. A metastable state 
could also exists for TV > iVcr, but estimating its lifetime 
is a more demanding calculation. One should also explore 
the energy barrier associated with other possible decay 
processes. Further work is planned in this direction. 
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FIG. 1. He density profiles (solid fines) in the z — 
plane of drops with N= 40, 100, 200, 300, 400, 500, and 1000 
having a vortex fine along tfie z-axis. The density profiles of 
vortex- free droplets are also shown (dashed fines). 



FIG. 2. The vortex energy AEv{N) (dots) in pure ''Hejv 
drops. The line is a fit obtained using Eq. (^) . Open symbols 
are the vortex energies AEy{N) in doped droplets. 



FIG. 3. Density distribution in the xz-pl&ne for a 
HCN-He5oo droplet hosting a vortex along the z-axis. Lengths 
are in A. Contour lines are drawn bewteen 10 equally spaced 
intervals of density, where white is for density less than 
5 X 10"^ A'^ and black for density higher than 5 x 10~^ A~^. 
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FIG. 4. Solvation energy Sx{N) for A= Ne, Xe, SFg, and 
HCN. The lines have been drawn to guide the eye. 

FIG. 5. Excess energies 5'x+v(A) (dots) for Ne, Xe, SFe, 
and HCN. The triangles represent 5'x (A) + A£v(A). Lines 
have been drawn to guide the eye. 
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